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ABSTRACT 

We present SHARC-2 350 /im data on 20 luminous z ^ 2 starbursts with 5i.2mm>2mJy from the Spitzer- 
selected samples of Lonsdale et al. and Fiolet et al. All the sources were detected, with 5350;^m>25 mJy for 18 
of them. With the data, we determine precise dust temperatures and luminosities for these galaxies using both 
single-temperature fits and models with power-law mass-temperature distributions. We derive appropriate 
formulae to use when optical depths are non-negligible. Our models provide an excellent fit to the 6 fim- 
2 mm measurements of local starbursts. We find characteristic single-component temperatures Ti ~35.5±2.2 K 
and integrated infrared (IR) luminosities around 10'^^^" 'Lq for the SWIRE-selected sources. Molecular 
gas masses are estimated at ~4x 10'" Mq, assuming K850pm=0.15 m^kg"' and a submillimeter-selected galaxy 
(SMG)-like gas-to-dust mass ratio. The best-fit models imply >2kpc emission scales. We also note a tight 
correlation between rest-frame 1 .4 GHz radio and IR luminosities confirming star formation as the predominant 
power source. The far-IR properties of our sample are indistinguishable from the purely submillimeter-selected 
populations from current surveys. We therefore conclude that our original selection criteria, based on mid-IR 
colors and 24 /im flux densities, provides an effective means for the study of SMGs at z~L5-2.5. 
Subject headings: galaxies: high-redshift — galaxies: ISM — galaxies: photometry — galaxies: starburst — 
infrared: galaxies — submillimeter: galaxies 
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1. INTRODUCTION 

Star-forming galaxies release a significant fraction of their 
energy output at infrared (IR) and submillimeter wavelengths 
(A^5-1000/im in the rest frame). The luminosities of star- 
bursts, with star formation rates over 100 M© yr"', are almost 
exclusively carried at IR wavelengths. Whereas luminous and 
ultra-luminous infrared galaxies (LIRGs and ULIRGs) are ex- 
tremely rare ir i the local universe, th e y are prevalent at h igher 
redshifts (e.g. ' Le Floc'h et ai]l2005t ICaputi et al.ll2007 ). By 
studying this population we can learn about the star formation 
history of the universe, and because the starbursting is trig- 
gered by merger events, we can also probe models of structure 
formation and halo dynamics through the ages. 

Around 1 mm wavelengths it is possible to find similar star- 
bursts across much of the volume of the universe, because 
the dimming of radiation from increasing distances is coun- 
tered by the steeply rising energy spectrum between 2 mm and 
200 /im in the rest frame. Thus, an IR galaxy would produce 
nearly the same (sub)millimeter flux regardless of its distance 
in the range of z 0.5-10. 

The nearly bias-free selection, together with the fundamen- 
tal d esire for studying starbursts, rnake submilhmeter sur- 
veys (ICoppin et al.l2006l:IWeiss et al.l2009l : [Austermann et al.l 
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'2010'; 'Iviso n et"ani201 0') especially relevant. However, while 
finding submillimeter-selected galaxies (SMGs) is relatively 
straightforward, studying them in any detail has proved diffi- 
cult. This is due to two factors: they are faint at other wave- 
lengths, and the re are often multiple possible o ptical/NIR 
counterparts (see lPope et alj2006tlYounger et al.l2009i) due to 
the poor spatial resolution of most (sub)millimeter telescopes 
(typicafly 10"-30" FWHM). 

As a result, most spectroscopic redshifts have been ob- 
tained for SMGs with radio-detected coun terparts provid- 
ing the requ ired positional acc u racy. Thus, IChapman et alj 
( 120031 12005h and iKovacs et all (120061) discovered that most 
SMGs lie around a median redshift of 2.3, are extraordinar- 
ily luminous (lO'^-lO'-' L©) and have large molecular gas 
reservoirs (10'"-10" Mq). The close correlation between IR 
and radio luminosities of SMGs resembles that of local star- 
forming galax ies (Helo u et al.1 119851: ICondon & BrodericH 
119911: ICondonI ri992t lYun et al.1 l20'0lh . implying that star- 
formation is the principal power source. 

Despite these successes however, the number of SMGs with 
spectroscopic redshifts is only around a hundred, and less 
than half of these are characterized in the far-infrared (FIR). 
Furthermore, the radio-undetected SMG population is almost 
completely unexplored. A truly unbiased understanding of 
SMGs requires improved selection methods. 
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Table 1 




Sample Selection Criteria 


B and/Instrument 


Criterion 


IRAC 




MIPS 24 


S24^.m>400AtJy 


MAMBO 


Si. 2mm > 2mJy 



Note. — Summary of the sel ectio n criteria for the MAMBO-detected 
5.8 fim peaker samples of iLOf^ and lF09t 



1.1. SMGs and Spitzer 

A major step forward in understanding SMGs was provided 
by Spitzer. With its sensitive mid-IR imaging and spectro- 
scopic capabilities, Spitzer detected most SMGs in multiple 
mid-lR bands, which can provide accurate photometric red- 
shifts dz/(l +z)~0.1 (Pope_el_al. 2006, Lonsdale et al. 2009). 
Furthermore, many Spitzer-selected objects are predic ted t o 
be submillimeter bright (Lonsdale et al. 2009; hereafter'L09'). 
Unfortunately, the identification of mid-lR counterparts to 
SMGs is plagued by problems similar to those of optical as- 
sociation (Pope et al. 2006), and so the reliance on radio as- 
sociations remains. 

Selecting strongly luminous starbursts at z^2 from Spitzer 
data can overcome these problems, and the resulting samples 
are expected to contain a large proportion of SMGs within 
them. Thus, a two-step process has been suggested. First, 
selecting sources with peak Spitzer-lnf mied Array Cam- 
era (IRAC) flux densities in the 5.8 /im channel should re- 
move galaxi es with a mid-IR luminous active galac tic nu- 
cleus (AGN; IWeedman et aLll2006l: iFarrah et aLllfool . leav- 
ing objects with a clear rest-frame 1 .6 /xm opacity minimum 
at zw2. The required IRAC detections ensure that stellar 
masses are large, especially wh en applied to limited sensi- 
tivity samples such as SWIRE (iLonsdale et al. 2003). Sec- 
ond, a bright 24 /im flux density cut (S24fj,m>^0Q /iJy) favors 
starbursts at the same redshift as the strong rest-frame 7.7 fim 
polycyclic aromatic hydrocarbon (PAH) emission feature is 
redshifted into the band. Redshifts from InfraRed Spectro- 
graph (IRS) spectra or from IRAC photome try confirm that 
most sources selected in this way lie at z—2 (IWeedman et al.l 
I2OO6I N. Fiolet et al., in preparation). The eff ective ness of 
such a two-pronged selecti on w as confirmed bylLQpl and Fi- 
olet et al. (2009, hereafter lF09h . who detected a significant 
fraction of these sources at 1.2 mm using the MAMBO cam- 
era at the IRAM 30m telescope (i.e., 5i.2|™,^2mJy). Clearly, 
the proposed Spitzer selection (see Table[lJi yields distant lu- 
minous starburst galaxies without radio preselection. 

A key question is to what extent the SWIRE-selected sam- 
ples are in fact representative of the purely submillimeter- 
selected population. A straightforward way to test this is to 
constrain where the peak of the IR emission lies, thus bridg- 
ing the gap between the cold dust detected by MAMBO, and 
the hot dust and PAHs detected by Spitzer. Such constraints 
can be used to derive effective dust temperatures and accu- 
rate IR luminosities, for comparison with the overall SMG 
population. Because the sample consists of galaxies at z^2, 
and dust temperatures are expected in the range of 30-40 K, 
the 350 /im band offers an ideal opportunity for providing the 
constraints we seek. 

For our study, we targeted 12 galaxies from |F09| and an- 
other 8 from L09. The former is an almost complete sam- 
ple of sources from 0.5 deg^ satisfying the criteria of Table [T] 



Table 2 

Variable Calibrators 



Name 


Date 


"obs 


(■ly) 


Cross 
Calibrators 


CIT6 


2007 Jan 17 


1 


3.14±0.25 


Arp220 




2007 Apr 17 


2 


3.12±0.11 


Arp220, Vesta 




2008 Jan 18 


5 


2.21 ±0.08 


Arp220 




2009 Jan 22-25 


16 


2.93 ±0.04 


Arp 220, Ceres 


3C273 


2007 Apr 19 


1 


4.08 ±0.30 


CIT6 




2008 Feb 23 


2 


3.82±0.21 


Arp 220 




2008 Feb 24 


2 


3.84 ±0.20 


Arp 220 




2009 Jan 23 


1 


5.04±0.37 


Arp 220 




2009 Jan 24 


3 


5.46±0.22 


Arp 220, Ceres 




2009 Jan 25 


2 


5.25 ±0.25 


Arp 220, Ceres 



Note. — The columns are: the name of the variable calibrator source, the 
period of measurements, the number of measurement, the derived 350 /im 
flux densities, and the stable cross-calibrator sources used. 



whereas the|L09l sources are drawn arbitrarily from ~30deg^. 
To enhance our chances of detection at 350 /im, we selected 
objects with clear MAMBO detections, with an additional 
preference toward the highest 1 .2 mm flux density measure- 
ments. 

The observations are described in Section |2] and the results 
are summarized in Section [3] where we also discuss one of 
the submillimeter sources, which appears to be a close asso- 
ciation of several IR galaxies. In Section |4] we develop ap- 
propriate spectral energy distribution (SED) models to inter- 
pret the measurements. With these models, we characterize 
the sample in Section |5] deriving precise dust temperatures, 
dust/gas masses and luminosities, and with them we verify 
the (far-)infrared to radio correlation. 

We assume a cosmology with //o=71 kms"' Mpc"', 
rj/ M=0.27 andJ7A=0.73, and calculate distances as prescribed 
bv lira([T999h . 

2. SHARC-2 350 /im OBSERVATIONS 

350 fim observations were carried out at the 10.4m Caltech 
Submillimeter Observatory (CSO), using the SHARC-2 Cam- 
era (iDowell et al . 2003). Data were collected during five suc- 
cessful observing runs between 2007 January and 2009 Febru- 
ary. The typical weather conditions during observing ranged 
from good to marginally useful, with in-band line-of-sight 
opacities t~1-2, corresponding to precipitable water vapor 
levels around 0.7-1 .3 mm. 

We took advantage of the op en-loop actuated D ish Surface 
Optimization System (DSOS; Leonget al.ll2006l) of the tele- 
scope to maintain the optimal beam qualities at all elevations. 
As a result, we expect our telescope beam efficiencies to be 
relatively high and limited mainly by the 12 /im surface rms 
of the CSO panels. 

Our observing strategy was to modulate the telescope pri- 
mary with small, field -of-view-sized open Lissajous patterns 
(lKovacsll2008Rl2006h . with average scanning speeds around 
a beam per second. These patterns facilitate the separation 
of the astronomical source from the various correlated sig- 
nals (of the atmosphere or the instrument). At the same time, 
the pattern provides uniform coverage over the field of view, 
the source remains exposed on the array at all times, and the 
scanning mode does not require extreme maneuverability of 
the telescope structure. 

In total, we spent 37.5 hr on the 20 objects observed. We 
aimed for 5 a detections, or limiting rms noise levels of 10- 
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Table 3 

Summary of the FIR and Radio Observations 









Offset 




5(24 /im) 


5(350 /xm) 


5(1.2 mm) 


5(20 cm) 


5(50 cm) 


5(90 cm) 


ID 




lAU Name 


(arcsec) 


z 


(MJy) 


(mjy) 


(mJy) 


(AtJy) 


(AtJy) 


(AtJy) 


LH-02 


SWIRE. 


_J103639. 57+575346.6 


+2.2, 


+1.0 


1.93" 


838±15 


54.8±5.8 


5.29 ±0.85 








LH-06 


SWIRE 


J103837.03+582214.7 


+3.2, 


-0.0 


1.68* 


1072 ±19 


45.5±7.1 


3.83±0.84 








LH-03 


SWIRE 


J104313. 33+574621.0 


+1.5, 


-0.2 


2.67" 


722 ±17 


62.9 ±10.2 


3.79 ±0.76 








L-1 


SWIRE 


J104351. 16+590057.9 


+2.7] 


-1.0 


2.26" 


722 ±17 


34.1 ±6.7 


2.95 ±0.66 


77 ±9.0 


142 ±19 


327 ±74 


L-9 


SWIRE 


J104440.25+585928.4 


+0.6, 


+0.4 


2.01'' 


674±19 


57.0±15.1 


4.00 ±0.55 


117±9.2 


195 ±22 


266 ±69 


L-11 


SWIRE. 


_J104556.90+585318.9 


+1.3, 


+0.8 


1.95' 


650±18 


39.7±5.9 


3.08 ±0.58 


315±4.1 


425 ±43 


651 ±72 


L-I4 


SWIRE. 


_J104638. 67+585612. 6 


+1.6, 


-0.1 


2.07' 


611±17 


31.9±4.9 


2.13±0.71 


160 ±6.0 


322 ±33 


426 ±70 


L-15 


SWIRE. 


J 1 04656.47+590235 .5 


+6.3, 


-0.9 


1.89' 


419± 18 


26.9±4.2 


2.36 ±0.62 


69 ±3.7 


181 ±21 


258 ±69 


L-I7 


S WiKb 


J 1U4 /04.9 /+5y2jM.j 


-1.2, 


+0.4 


1.99' 


647 ±18 


39.6±8.5 


2.24 ±0.64 


341 ± 33 


756 ±78 


1 rid 1 "TC 

1052 ± 75 


L-20 


SWIRE. 


_J104717.96+590231.8 


+4.7, 


-1.9 


2.10' 


617±19 


49.7±6.5 


2.66 ±0.78 


51 ±4.7 


166±19 


264 ±72 


L-21 


SWIRE. 


_J104718.63+584318.1 


+0.3, 


+5.3 


1.78'' 


447±16 


41.4±7.4 


3.09±0.81 


138±36 


176±39 


350±76 


L-22 


SWIRE. 


_J104720.49+591043.6 


+2.5, 


+1.5 


2.57'' 


435 ±16 


18.9±5.7 


3.41 ±0.73 


101±7.1 


238 ±26 


447 ±70 


L-23 


SWIRE 


J104726.44+585213.3 


+4, 


-3 


2.38'' 


434 ±17 


19 ±9' 


3.13±0.86 


86±7.2 


172±21 


242±71 


L-25 


SWIRE 


J104738.32+591010.0 


+1.2, 


-0.3 


1.96' 


723 ±17 


31.8±5.9 


2.56 ±0.74 


69 ±9.0 


166 ±20 


348 ±72 


L-27 


SWIRE 


J104744.59+591413.4 


+0.9, 


-2.4 


2.20' 


523 ±18 


28.7±6.5 


2.48 ±0.73 


77 ±14 


106±18 


<219 


LH-01 


SWIRE 


J105007.26+571651.0 


+1.8, 


+2.3 


1.35" 


756±12 


46.6 ±7.6 


5.65 ±0.74 








ENl-OI 


SWIRE 


J160343.08+55 1735.8 


+2.9, 


+0.1 


2.37" 


699 ± 14 


42.7 ±8.2 


2.81 ±0.78 








EN 1-02 


SWIRE. 


_J160440.45+543 103.0 


-0.5, 


-7.3 


1.93" 


772 ±18 


44.3 ±9.1 


2.46±0.81 








EN 1-04^ 


SWIRE 


J161658.14+535319.3 


+0.2, 


+0.6 


2.00" 


657 ±15 


27.6 ±10.2 


2.82 ±0.99 








EN2-01 


SWIRE. 


_J163734.44+415151.4 


+0.9, 


+3.8 


1.85' 


8 14 ±20 


30.4 ±5.7 


2.52 ±0.70 









Note. — Source IDs prefixed with LH or EN are from'L09', whereas those lettered with L are from'F09'. SHARC-2 detection offsets are with respect to the 
SWIRE catalog coordinates. 1.2mm flux densities are from L09 and F09, 20cm radio data from Owen & Morrison (2008), 50cm value from F. Owen et al. (in 
preparation), while 90 cm points are from cOwen et al.ii20Q9.1 . The quoted rms uncertainties do not include the uncertainties in the calibration. 

" Photometric redshift from LOSj. 

IRS Spectroscopic redshift from lL09l 
' IRS Spectroscopic redshift from N. Fiolet et al. (in preparation). 
"' Photometric redshift from F09. 

' The 350 flux density of L-23 is from the undeconvolved image close to the instrument resolution (middle panel of Figure[T). 

f SWIRE_J161658. 14+535319.3 is the only source in our sample with a MIPS 160 fj,m detection of 81±4mJy. However, this datum is not included in our fits. 



12mJybeam"'. Pointing and flux density calibration mea- 
surements were taken frequently, at least once per hour during 
the observing period. Our principal flux density calibrators 
were Mars, Arp220 (9.93±0.11Jy, based on our own cal- 
ibration against planets and asteroids), and the minor plan- 
ets Ceres, Juno and Vesta, for which appropriate temperature 
models have been derived by C. D. Dowell and are available 
through the SHARC-2 Web siteQ In addition, we relied on the 
evolved star CIT6, and the quasar 3C273 as additional cali- 
brators, owing to their relative proximity to the targeted fields 
in the Lockman Hole. To counter the long-term variability of 
CIT6 we measured it against our main calibrator sources dur- 
ing each observing run, while 3C 273 we cross-calibrated on a 
daily basis (see Table|2]l. Based on the residual scatter of cali- 
brator flux densities, we expect a 350 /.tm calibration accuracy 
around 8% for our science targets. 

We checked and adjusted the focus settings several times 
during each night, more frequently near sunrise and sunset 
when variations tend to be greatest. To avoid the unnecessary 
smearing of our sources due to systematic pointing errors, the 
pointing data of each run were carefully analyzed for trends. 
Corrections were applied as required at the time of reduction. 
The final pointing accuracy of the SHARC-2 data is expected 
at ~3" rms. 

We r educed the data using the CRUSH software (1.99-a5; 
iKova cs 2008a, 2006), which estimates and removes the cor- 
related atmospheric and instrumental signals, solves for the 
relative detector gains, and determines the appropriate noise 
weighting of the time streams in an iterated pipeline scheme. 
To maintain calibration accuracy, the inherent filtering of 

* |www ■ submm . caltech . edu/~sharcl 



point-source signals by the reduction steps is carefully esti- 
mated and corrections are applied internally by CRUSH. 

3. RESULTS 

We detected all 20 of the targeted sources. In a few cases, 
the detections are tentative. The measured 350 /im flux densi- 
ties are listed in Table [3] along with all prior measurements 
relevant to this work. We reached rms sensitivities of 4- 
lOmJybeam"' in all but one field. 80% of the 350 /im de- 
tections are robust (>4 a), two are 3^ a and another two are 
marginal (<3(t). However, even the lowest signal-to-noise 
detection is significant (with a confidence level >95%) given 
the small, beam-sized search areas surrounding the IRAC (and 
nearby radio) positions. For the most part, the SHARC-2 po- 
sitions are in excellent agreement with those from IRAC, with 
deviations at <3" rms. 

3.1. Submillimeter Multiplets 

There is, however, the curious case of L-23 (see Figure [T]), 
with an elongated 5.7 cr SHARC-2 peak 13".5 from the IRAC 
and radio positions. The discrepancy is too large (>4(t) to 
be a statistical pointing error Coincidentally, we have a large 
number of pointing measurements both before and after ob- 
serving the galaxy in question, all reinforcing an absolute 
350 /im astrometric accuracy better than 3" rms. Thus, the 
single SHARC-2 peak is inconsistent with a source at the 
IRAC and radio positions. 

A substructure of the peak is revealed if we smooth less 
than the optimal amount for point source extraction. We find a 
banana-shaped feature dominating the 350 /im emission to the 
east, and a fainter ~2 a compact peak close to the expected 
position. The position and point-like nature of the fainter 
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Figure 1. 45" x 45" postage stamps around L-23. Left: the beam-smoothed image used for source extraction at ~ 11" resolution shows a single peak east of the 
expected SWIRE source position (dashed circle indicating typical search area). Center: near the diffraction-limited resolution of SHARC-2, the peak resolves 
into a banana-shaped feature and a hint of a point source close to the IRAC position. The image has been slightly realigned to match the elongated feature with 
two nearby IRAC galaxies (yellow squares). The nearby radio sources are shown with red diamonds. The search area is reduced as a result of the reahgnment. 
Right: a deconvolution suggests up to four distinct components (lettered A through D), all within the typical beam areas of current submillimeter surveys. Three 
of the peaks align well with IRAC sources in the field. Contours are -15,0, 15, 30, and 45 mjy beam"' for the panels on the left and center, and -7.5, 0, 7.5, 15, 
and 22.5 mJy beam"' for the image on the right. 



Table 4 

350 fim Decomposition of L-23 



ID ^2000 "52000 S35o^m Association 



A 10''47"'28M -h58°52'20" 25mjy IRAC/M1PS24, radio 

B 10'"47™27^6 -l-58°52'24" 17mjy Unknown 

C 10''47"'28M -l-58°52'15" 14mJy IRAC/M1PS24, radio 

D 10''47"26'.9 +58°52'14" lOmJy L-23 



fact that it is a triplet or a quadruplet suggests that unresolved 
doublet, or even triplet (sub)millimeter sources may abound. 
The ^20% fraction of t he radio-identified SMGs with mul- 
tiple radio counterparts (llvison et alj|2007l) could be an indi- 
cation of the frequency of unresol ved (sub)millim eter multi- 
plets, which is also consistent with I Ara vena et all ( 2010). As 
such, the ubiquitous presumption that (sub)millimeter sources 
are individual galaxies may not be entirely justified. 

If multiple luminous galaxies are often behind the low- 
resolution (sub)millimeter peaks, the implications would be 
wide ranging. Some or all of the IRAC sources com- 
monly de tected close to an uncertain (sub)millimeter posi- 
tion (e.g. iPope et alJ 12006: Younger etal. 2009) could each 
contribute comparably to the measured (sub)millimeter flux 
densities. A commonality of multiple ts would imply that 
(sub)millim eter source counts may be s teeper than what 
SHAD ES (iCoppin et al.1 l2006i) . AzTEC dAustermann et alj 
HnH), or the LABOCA deep-field survey (Weiss et al. 200% 
found; by counting multiple faint sources as a single bright 
one, these surveys would slightly underestimate the abun- 
dance of faint galaxies while substantially overestimating the 
number of bright ones (i.e., ULIRGs and hyper-LIRGs). Sub- 
millimeter surveys by Herschel (especially at 250/500 /im) 
would be similarly affected. A strong clustering of SMGs at 
separation angles <20", i.e., below the smallest scales that 
current (sub)millimeter surveys can resolve, would be im- 
plied. Unresolved multiplets could partially explain why as 
many as half of the SMGs in the literature have no identifi- 
able radio counterparts: as opposed to having one detectable 
counterpart, there may be a number of associated faint radio 
sources below the detection level. 

4. SPECTRAL ENERGY DISTRIBUTIONS 

The addition of the 350 jim data to the existing measure- 
ments provides a powerful c onstr aint on th e FIR SEDs of the 
SWIRE-selected samples of|LOland|FS It allows, for the 
first time, an accurate characterization of the dust tempera- 
tures and luminosities of individual sources in the sample, 
provided that physically realistic models of emission are fitted 
to the observations. 

Before the emergence of details on the Wien-side shape of 
SEDs (e.g., from Spitzer), it was common practice to rely 
on single-temperature dust models. Extra temperature com- 
ponents were added as necessary to fit data at wavelengths 
<100/zm . However, such models were clearly very naive. 



Note. — The most likely decomposition of the blended 350 fj,m emis- 
sion from the neighborhood of L-23 I F09), ordered in decreasing estimated 
brightness. The SWIRE/MAMBO source we identify as component D. The 
uncertainties in the deconvolved flux densities we estimate at around 5 mJy. 
Positions are estimated to be accurate to <2" rms. The positi ons o f the asso- 
ciated IRAC/MIPS and radio sources are discussed in Section IxT] 

component is consistent with both our pointing accuracy, and 
the expected size (<!") of a galaxy at z=2.38. Therefore, we 
identify this as the 350 /im counterpart to L-23. 

Deconvolution of the elongated feature suggests three 
distinct components. Two of these align well with 
nearby IRAC and MIPS 24 /im peaks (at 10H7'"28M4, 
-l-58°52'20".9 and 10H7'"28^04, -l-58°52' 14".l) and radio 
sources (at 10H7'"28^45, -l-58°52'20".9 and 10H7"27^97, 
+58°52' 14".l), while the third peak is perhaps a super- 
resolution artifact. The positions and approximate brightness 
of the deconvolved peaks are summarized in Table H) 

It is not unusual to find resolved pairs of submillime- 
ter source s with ~'20"-30" separation (Blainetal. 2004; 
iGreve etal.. 2004: Farrah et al. 2006; Weiss et al. 2009). Such 
pairs are relatively infrequent and correspond to massive dark 
matter halos (>10'^ M©). However, the clustering of SMGs 
on smaller scales is not well constrained. Therefore, the 
blending of as many as three or four galaxies into a single 
350 /im detection peak (with total 5'35o^m~60mJy) is inter- 
esting. 

We often assume that each deep-field (sub)millimeter 
detection corresponds to just one luminous galaxy, with 
unique counterparts at optical, IR, or radio wavelength. 
Granted that the clustering of SMGs on scales <15" is 
not well known, it seems reasonable given the generally 
low density ('^400/deg^) of sources in the current, shallow 
(sub)millimeter surveys. However, our multiplet puts this as- 
sumption to the test: the 350 /im components that we identify 
around L-23 would constitute a single source under the spatial 
resolution (15"-30") of current (sub)millimeter surveys. The 
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Among the many m ore physically motivated trea t ments is the 
pioneering work by iDesert. Boulanger & Pugetl (11990b who 
considered dusty environments exposed to a range of radiation 
fields. Their framework has been fine-tuned to better match 
observations, especially of the aromatic features in the mid- 
IR (e.g. Dale et al. 2001; D ale & Helou 2002). These models, 
fitted to well-studied local galaxies, serve as typical templates 
for interpreting Spitzer data. 

Our goal is to derive accurate and robust FIR characteri- 
zations for our sample of galaxies. In so doing, we wish to 
move beyond a hodge-podge of temperature components but 
we also do not need the full complexity of the models with 
mid-IR detail. Thus, we fitted both single-temperature mod- 
els (which offer direct comparisons to many previous stud- 
ies) and power-law temperature distributions for a more accu- 
rate description. In the following sections, we shall discuss 
the models and the formulae we use for deriving the essential 
quantities. 

4. 1 . Single-temperature Dust Models 

The simplest model of the thermal IR emission assumes 
that it can be described as a singular gray-body at a physi- 
cal temperature T in the rest frame. Such single-T flux den- 
sities S\ are best expressed in terms of the observables: the 
measurement frequency j^obs, and observed-frame temperature 
robs = r/(l+z), as 



(1) 



The other terms are the physical solid angle dVl, the opti- 
cal depth r, Planck's blackbody law B^, and magnification m 
(for gravitationally lensed objects). The physical solid angle 
is related to the projected source area A by the angular-size 
distance Da as dVl =A/D\. The optical depth can be esti- 
mated as the ratio of the total emitting cross-section of dust 
particles to the projected area of emission. As such, optical 
depths must be evaluated in the rest frame of the emission. 
Adopting a typical power-law frequency dependence for the 
emissivity for dust particles, optical depths can be expressed 
in the rest frame of the galaxy as 



1/3 



_Md_ 

dnDl' 



(2) 



Here, J^r = (1 +z)t'obs is the rest-frame frequency, is the 
total dust mass at an angular-size distance of Da, to is the 
characteristic photon cross-section to mass ratio of particles at 
frequency i^q. The spectral emissivity index /3 can be related 
to the fracta l dimension D of the dust particles as (3 = D-l 
(lYang|2007h . Thus, /3 values are expected in the range of 1-2, 
depending on the shape of particles. There have been several 
attempts for determining appropriate k values for intersteflar 
dust particles (e.g. Hildebrand 1983: Krugel. Steppe & ChiTiil 
[7990: Sodroski etal. J997- James et al. 2002t ). The resulting 
scatter of values reflects the difficulty of making such mea- 
surements. In this work, we assume Kg50um=0.15m-kg"', 
which matches the value used by Kovacs et al. (2006) and is 
based on an average 125 fim value (Dunne, Bales & Edmundsl 
|2003) assuming /3=1.5. As the strength of emission depends 
on the Ki,Md product, the inherent uncertainties in k ulti- 
mately affect the accuracy to which dust masses can be de- 
termined. At the same time, the temperatures and luminosi- 
ties can be determined accurately and independently from the 
assumptions on the absorption efficiency. 



We can calculate single-temperature IR luminosities Li by 
integrating the de-magnified observed flux density at luminos- 
ity distance Dl=Da(1 +z)^ from the source: 



LiiT)=4TT Dim-'- 



Siiv, T) dv 



The integration can be performed analytically in the opti- 
cally thin limit (t<^;1), giving 



Li(r) = r(4 + /3)C(4 + /3) 



Stt/i 



kT 



Aeff, 



(3) 



in terms of the effective surface area of emission Aeff, 
Planck's constant h, the Boltzmann constant k, the speed of 
light c, and Riemann's F and ( functions. The same expres- 
sion is also applicable for the case of optically thick emission 
(r^l) by setting /3=Q. The effective area of emission Aeff, in 
the optically thin and thick limits respectively, is given by 



MdKQikT /hvQ)l^ ifr<l 



ieff ' 



dQDl 



if T > 1 or /3 ^ 



(4) 



Given the relation of the projected surface area to total emit- 
ting area (7rr^ versus A-nr for a sphere), our expression for 
the optically thick limit recovers the Stefan-Boltzmann law 
of blackbody luminosities. In the optically thin limit, Aeff is 
the physical emitting area at the frequency v,—kT /h. From a 
computational point of view, it is practical to approximate the 
Riemann F and C, functions in Equation (O for the range of 
typical (3 values: 



F(4 + /3)C(4 + ^) w 6 .45 exp( 1 .24/3 + 0. 1 1 /S^) , 



(5) 



which is accurate to about a percent in the range /3~0-2, 
Interstellar dust emission is mostly optically thin at FIR 
wavelengths (A>100/im). Only the more extreme galaxies 
and quasars, or strongly lensed objects, concentrate dust to 
such an extent that optical depths can become significant over 
the bulk of the emission spectrum. This is a concern, given the 
extreme nature of our objects. We therefore would like to also 
consider the effect of intermediate optical depths in deriving 
accurate luminosity measures. 

4.1.1. Intermediate Optical Depths 

We can derive approximations for the integrated luminos- 
ity for the general case of non-negligible optical depths. Our 
approach is based on the observation that most of the lumi- 
nosity is radiated at frequencies close to the gray-body SED 
peak (see Equation ([T])), which is observed near 



i^pk«(3 + /3) 



kT, 



obs 



(6) 



Around the peak, we can approximate the frequency de- 
pendence of the emissivity, i.e. e = (1 -e~^), with a power law 
e ~ epk(i'/t'pk)'^°", where, by definition, an effective spectral 
index /3eff can be calculated as 



'eff = 



91oge 



V de 
e dv 



(7) 



d\ogv 

Substituting the optical depth from Equation (|2]) and evalu 
ating at the peak Tp^, we arrive at 



/3efr = /3 



(8) 
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Table 5 

FIR and Radio Properties of Local Starbursts 

















logMrf 


logL 








Name 


(Mpc) 


Mr 




lx,-l 


(K) 


7 


(Mq) 






q 


a 


Arp220 


77.6 


15 


7 


0.92 


36.4±1.2 


7.7(2) 


7.92(3) 


12.16(2) 


0.68 


2.47(4) 


0.25(4) 


IRAS 20551-4250 


188 


10 


1 


0.99 


37.9±2.4 


7.2(3) 


7.33(9) 


12.09(3) 


0.19 


2.61(9) 




IRAS 22491-1808 


348 


14 


1 


0.68 


39.0±2.0 


7.6(3) 


7.62(8) 


12.20(3) 


0.38 


2.91(7) 




M82 


3.57 


11 


9 


0.99 


35.4±2.0 


7.1(3) 


6.12(7) 


10.86(2) 


0.12 


2.39(8) 


0.65(5) 


NGC 6090 


127 


11 


1 


1.81 


28.0±1.9 


6.9(3) 


7.44(8) 


11.67(2) 


0.15 


2.33(12) 




NGC 6240 


105 


21 


7 


1.07 


32.3±1.1 


7.1(1) 


7.56(4) 


11.96(2) 


0.25 


1.84(4) 


0.85(3) 



Note. — The FIR and radio properties of local starbursts, from the multi-T fits to a selection of available IR continuum data in the 6^m-2mm range, with 
13=1.5 and assuming typical emission diameter of 1 kpc (or 300 pc for M 82). The synchrotron component is modeled as a simple power-law SiadO<:'^ Columns 
are: the object name, luminosity distance Dl, the number of IR and radio flux density data (Mr and Mad) used, the residual scatter \xr\ around the fit, the 
temperature Tc of the dominant cold component, the mass-temperature index 7, dust masses Mj, integrated IR luminosities L, the optical depth r around the IR 
emission peak u^^^^, the radio-(F)lR congelation constant q, and synchrotron spectral index a. Uncertainties are 1 cr total eiTors. These are indicated in brackets 
for 7, logMrf, logL, q, and a, and refer to the last digits in the quoted values. Dust masses assume K(850/im)=0. 15 m^ kg"' . Plots of these models and the data 
used for deriving them are available in the online edition. 



Thus, the frequency dependence of the emissivity near the 
peak of the emission can be characterized by an effective 
gray-body index /3eff ■ This means that we can rely on Equa- 
tion (O to provide an approximation to the luminosity func- 
tion, simply by replacing f3 with /3eff above together with an 
effective emission area: 



leff ' 



(i-g-^p") dnol 



(9) 



(3-i-/3)ft'f 

It is possible to refine this approximation by seeking si- 
multaneous solutions to Equations (|6]i and dHJ with /3— )>/3eff 
in Equations ^ and 

4.2. Power-law Temperature Distributions 

Single-temperature models can be useful, especially for of- 
fering comparison to earlier studies that analyzed data in this 
way. However, it is clear that such models fall short on the 
Wien side of the thermal emission spectrum (dotted line in 
left panel of Figure |2]i. Measurements at wavelengths shorter 
than the peak (typically ~200 /im) by IRAS, ISO, or Spitzer 
all indicate that flux densities fall far less dramatically than 
single-T gray-body models imply. Temperature estimates are 
biased by the inclusion of data on the Wien side (at <200 /im 
in the rest frame), and luminosities are systematically under- 
estimated. The errors can be significant depending on the ac- 
tual spectral profile, and the data used for fitting. This is a 
concern, since the SHARC-2 350 /im point translates to rest 
frame wavelengths of 90-175 /im at redshifts z^l-3. 

For improving our models, we look for clues in the 
same local starbursts that yield t he oft-used templates for 
the in terpretation of Spitzer data (lBertall2005t iPoUetta et al.l 
120071) . Based on these galaxies we find that a power-law 
distribution of temperature components (dMd/dT oc T~^), 
with a low-temperature cutoff at T^, are in excellent agree- 
ment with measurements (Figure |2]i. At the same time, 
power-law distributions are easily justi fied in a physical con- 
text. The radiative-tr ansfer models of iDale et al.l (1200 Ih and 
iDale & Heloul (120021) are similarly conditioned on dust being 
distributed in a range of radiation fields U as dM/dU oc f/"". 
Their power-law index a is also the Wien-side spectral index 
and depends on the geometric configuration of dust clouds 
relative to the heating sources. Under the assumption of ra- 
diative dust heating the two models are approximately equiv- 
alent with 7 « 4+a+fieff. Therefore, we expect 7w6.5-7.5 for 
sources embedded in a diffuse medium (e.g., star formation), 
and 7^4-5 if absorption happens inside a dense medium at 



a well-defined distance from the source (e.g., for a molecular 
ring surrounding an AGN). 

We argue that a temperature distribution (as opposed to a 
radiation-field distribution) offers a more general treatment, 
because it can also account for non-radiative heating of dust 
particles, e.g., by shocks or gravitational collapse. One needs 
only to remember that one of the critical functions of dust in 
star formation is to act as a coolant during the collapse of star- 
forming gas clouds. Thus, dust heating begins even before the 
UV-bright OB stars illuminate the interstellar medium. Tem- 
perature distributions can further arise from a distribution of 
individual molecular clouds in a galaxy, each with its own size 
and environment. 

Apart from the physical plausibility of power laws, a prac- 
tical advantage is that they add just a single variable 7 to the 
parameterization of single-T models. This is welcome, es- 
pecially since the number of SED measurements across the 
FIR band are generally few, even for many well-studied local 
galaxies, allowing robust constraining only for a handful of 
parameters. 

We can write appropriate expressions for predicted flux 
densities in terms of the single- T flux (Equation ([T]i) as 



5i(j^obs,7Lbs) T-'dT, 



(10) 



where once again Tobs = T /(l+z). The integration can be 
performed numerically with a number of discrete tempera- 
ture compone nts, or can be written in a closed form (see 
iKolbidfTOTOh using the incomplete Riemann zeta function 
Zi(7-1, hvr/kTc) and r(7-l). Luminosities can be ex- 
pressed in a closed form and are given by 



7-1 



(11) 



7-(5+/?eff) 

based on the single-T value L\ (Equation (O). Finiteness of 
luminosities requires that 7 be larger than 5+(i^ff. 

Our measure of the luminosity accounts for the full power 
radiated through thermal dust emission, in contrast to a more 
traditionally defined FIR or IR luminosities, which measure 
output inside a strictly defined wavelength ranges (e.g., be- 
tween 40 /im and 1000 /im). As such, our expression for the 
integrated luminosity can be applied to a wider range of envi- 
ronments (e.g., hot accretion disks) than definitions based on 
strict wavelength ranges. Additionally, L provides a fair proxy 
for the total bolometric luminosity of warm (>30 K) LIRGs, 
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^ ^cpeak ^ ^ ^cpeak 

Figure 2. Composite radio and IR continuum SEDs of local starbursts (left) and the SWIRE/MAMBO sample (right). The gray lines represent the best-fit gray- 
body models with power-law mass-temperature distributions using_7=7.2 for the local starbursts and 7=6.7 for the SWIRE-selected sample. Frequencies are 
normalized to the peak of the cold-component emission (Equation ^). A single-T model is shown in the left panel for comparison (dashed line), demonstrating 
the shortcoming of such models on the Wien side. The continuum contribution to the 24 fim flux density is estimated based on the average IRS spectra of eight 
IL09 galaxies. The average local starburst model with 7=7.2 is also shown in the right panel for comparison (dashed Hne). We also show the stacked average 
70 fim and 160 fim flux densities of F09 sources (black hollow squares with error bars) for a Active source with median properties (7'r=33.5 K, Mrf=7 X 10^ Mq 
at z=2.05).The sources with insufficient radio data assume Is = 0.75 for the radio spectral index. Individual SED plots for the local starbursts (part b) and the 
SWIRE-selected sample (part c) are available in the onHne edition. For part c, the 24 lijm data points reflect the expected continuum contribution to the in-band 
flux-densities, according to Figure 3, with an uncertainty of 0.25 dex. (An extended color version of this figure is available in the onhne journal.) 



where an overwhelming fraction of the total power output is 
absorbed and re-radiated by dust. 

4.2. 1 . Clues from Local Starbursts 

To check what /3 and 7 describe the local starbursts best, 
we fitted a selectior0 of published FIR continuum data un- 
contaminated by aromatic features or synchrotron emission, 
in the 6 /im-2 mm range, with power-law T distributions and 
assuming kpc-scale emission. For M82 we have assumed 
a characteristic emission diameter of 300 pc, in line with 
the observed size of its molecular ring. Our fitting routine 
performs a generic minimi zation based on the downhill 
simplex method described by iPress. Flannerv & Teukolskvl 
(|T986). The results are summarized in Table|5] 

Encouraged by the consistent 7 values within measurement 
error, we stipulate that the temperature distribution is possi- 
bly sufficiently universal for star-forming galaxies, so that we 
may proceed to fit 7 and /? globally, while keeping Tc and 
as individual parameters for each galaxy. We do this by nest- 
ing an inner minimization loop for Tc and Md for each ob- 
jects, in an external fitting of /3 and 7 for all. Thus we find that 
13 = 1 .53 ±0.04 and 7 = 7.22 ±0.09 describe the local starburst 
best, producing close to the expected level of residual scat- 
ter based on the measurement and calibration uncertainties. 
The resulting /3 is i n excellent agreeme nt with measurements 
both in the Galaxy ( Pupae et al.ll20()3l) and in the laboratory 
jAgladze et al.l l 1 996*). Under the assumption of radiative dust 
heating, the values of 7 and the characteristic opacities t>0.3 

^ We relied on the online NASA Extragalactic Database (NED) for flux 
densities. When more than one measurement was available for a given band 
we generally opted for the most recent one, or the one with the lowest mea- 
surement uncertainty. Additionally, we supplemented the data for M 82 in the 
350/xm-2mm range with Leeuw & Robson (2009) and Alton et al. 1 1999). 
For Arp 220, we also included our own 350 fim measurement of 9.83itO. 11 Jy 
in the analysis. 



on the Wien side imply radiation-field index a>2.35, in good 
agreement with the predic ted value of 2.5 for sources embed- 
ded in a diffuse medium CDale et al.ll200 1). The conclusions 
on (3 and 7 hold for a wide range of possible emission scales 
around the assumed values. 

4.3. Sample Characterization 

In characterizing our SWIRE-selected samples, we assume 
/3=1.5 based on our findings on local starbursts and work 
by others. For the size of the emitting region (i.e., for set- 
ting dflAD^), we assume an average diameter of 3kpc in 
our single-temperature fits. Such an extended starburst is un- 
precedented in the local universe, where the luminosities tend 
to originate from smaller kpc-scale molecular rings. Never- 
thel ess, the larger size s are reasonable for our sample (see, 
e.g. Farrah et alJl2008 ). where luminosities exceed the local 
examples by 1 or 2 orders of magnitude. The assumption 
on the extent of emission is critical only for interpreting the 
temperature parameter of the fit as a true measure of physical 
temperatures. The best-fit temperature values would rise if 
smaller emission diameters were assumed. At the same time, 
dust masses and integrated luminosities remain nearly inde- 
pendent of the size assumption, and thus offer robust charac- 
terizations. The resulting single-temperature fits to the sample 
are summarized in Table |6] 

The FIR properties derived from the single-T fits are based 
on the MAMBO 1.2 mm measurement and our SHARC-2 
350 /im data only. Shorter wavelength data (especially MIPS 
24 //m) were not included in these fits because these trace 
higher-temperature components and aromatic features, nei- 
ther of which are accounted for in the single- T dust model. 
We also do not include the 160/im datum of EN 1-04 for 
any in our modeling, because we find it inconsistent with the 
measurements at 24 /um and at 1 .2 mm. We suspect that the 
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uniquely detected 160 /im flux density of our sample may be 
contaminated by an unknown source close to EN 1-04. 

Even with all the caution we can summon, single- 
temperature models cannot yield accurate characterizations 
for reasons already disserted. We include these simple mod- 
els only to provide a comparison to other studies (especially 
iKovacs et al.l 120061) . which analyze FIR and (sub)millimeter 
data in this way. Fitting a distribution of dust temperatures to 
the data yields more realistic estimates of the dust masses and 
luminosities. Our best estimates of the underlying FIR prop- 
erties are summarized in Table |7] and shown in Figure |2] The 
temperature parameter Tc indicates the temperatures of the 
dominant coldest component, above whic h w armer compo- 
nents exist in lesser quantities (see Section l4!2T i. Therefore, it 
is not surprising that single-component fits, which measure an 
average temperature, should yield characteristic values above 
the cold cutoff temperatures Tc in a distribution. 

For our temperature-distribution fits, we have included 
the MIPS 24 /im data as guiding constraints on the short- 
wavelength tails of the dust continuum. It is clear, however, 
that the 24 jiva flux density is likely dominated by emission 
from aromatic features. However, warm dust and PAH likely 
cohabit the same interstellar environments. Thus, PAH emis- 
sion is almost certainly coupled to a comparably bright dust 
continuum. Consequently, one should anticipate close to con- 
stant PAH-to-continuum flux ratios for a class of star forma- 
tion dominated objects. Based on the stacked IRS spectra of 
IL09I we can estimate the typical continuum contribution to 
the MIPS 24 ^va measurement as a function of redshift (Fig- 
ure [3]i, with which we set appropriate 24 /im constraints for 
our fits. Recognizing, however that the continuum fraction 
may vary significantly from one galaxy to the next, we as- 
sign an uncertainty of 0.25 dex to these 24 /im guesstimates. 
Thus, we allow the 24 /im continuum fraction to span a full 
decade in brightness before the deviation becomes significant 




1 1.5 2 2.5 3 

Redshift 

Figure 3. Expected continuum contribution to the MIPS 24 /im measure- 
ments. Tlie top panel sliows tlie average IRS spectra of eight L09 sources 
(thin solid line), and as it would appear under the resolution of the MIPS 
24 /im band (thick solid curve). Based on the decomposition of this spectrum 
into starburst PAH features (NGC3079) and a power-law continuum (dashed 
line) contribution ~A- (see L09), we estimate the typical continuum contri- 
butions to the MIPS 24 /im flux as a function of redshift (bottom panel) for 
the SWIRE-selected samples of L09I and lF09 . which we use to derive appro- 
priate continuum constraints at 24 /im for our SED fits. 



(i. e., > 2g), while we assume that the average IRS spectrum 
of IL09I is representative of the sample as a whole. 

5. DISCUSSION 

Under the assumption of extended starbursts we find that 
our models deviate from the optically thin limit only mod- 
erately, with Tw0.7±0.2 near the emission peaks (under the 
assumed emission diameters of 3 kpc). This is good news for 
a straightforward comparison to the SMGs of iKovacs et al.l 
(2006), which were analyzed under the optically thin assump- 
tion. These optical depths are not much higher than the opac- 
ities of local starbursts on kpc scales. We may take this as a 
first indication that our guess on the extent of star formation 
is correct. 

The median 70 /xm and 160 /xm flux densities of 
2.56±0.38mJy and 14.63±4.05mJy, respectively, from the 
stacked images of the MAMBO-detected sources of F09 agree 
well with the multi-T models of our sources with median 
properties (i.e., r,.=33.5K, Mrf=10'*-^2MQ at z=2.05), within 
the expected measurement uncertainties (see Figure|2]l. These 
prove the general reliability of the 24 /im continuum contri- 
bution estimates. 

The use of the scaled 24 /im data further allows us to 
probe the range of 7 values and typical emission scales that 
are consistent with measurements. Accordingly, we find 
7=6.71 ±0.11 and average emission diameters around ~2kpc 
fit the data best globally. We can exclude diameters < 1.2 kpc 
with 95% confidence. Thus, we see our size assumption for 
the single-T fits sufficiently justified. At the same time, we 
adopt the most likely diameter of 2 kpc in our multi-T fits. 

The result on the sizes is significant because it provides 
an indirect measure of the true extent of the FIR emission 
in hyper-lumino us infrared galaxies (HLIRGs). Interferomet- 
ric imaging by Tacco ni et al.l (|2006) akeady hinted at similar 
scales based on slightly broadened CO emission profiles (ver- 
sus their synthesized beams). The broadening however can 
also result from distinct compact nuclei separated at the in- 
ferred distance (e.g., in a merging system), whereas the best- 
fit SEDs are only consistent with extended dust emission. The 
implied scal es, spanning kiloparsecs across, are in line with 
iFarrah et"ar ( 2008) and are a strong indication that dust heat- 
ing is distributed, which in turn supports star formation as the 
primary power source. 

F or star-forming galaxies, we expect 7«6.5-i-/3eff based 
on (iDale et al.ll2001h . With Tpk>1.2±0.4 near the emission 
peaks (corresponding to the most likely emission diameters 
of 2 kpc), it is not surprising that the SWIRE-selected sam- 
ple has 7 values below those of local starbursts. The rela- 
tively shallow temperature distributions in turn might explain 
the unusually high 24 /im to 1 .2 mm flux density ratios re- 
ported by L09 and F09. The implied radiation-field distribu- 
tion indices Q!>2.0-2.35 further support star formation as the 
primary heating source behind the extreme luminosities. 

The average dust mass in our sample is 7xlO'^Mo. All 
of our galaxies have masses near the mean value, within 
the typical 40% (i.e., 0.15 dex) measurement error, revealing 
only a hint of an underlying scatt er <20%. Assuming SMG- 
like gas-to-dust mass ratios of 54 ("Kovacs et al."2006), which 
are also typical of t he nuclear regions of nearby galaxies 
(ISeaquist et al.ll2004l) . we estimate an average molecular gas 
content for our sample in the neighborhood of 4 x 1 M© . 

Ou r a nalysi s reveals no statistical difference between the 
'L09' and IF09' subsamples. Table [8] offers a comparison of 
average properties between the subsamples, the combined 
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Table 6 

Single-temperature Characterizations 





T 


logMrf 


logLi 




ID 


(K) 


(Mq) 


(Lq) 




LH-02 


36.5 ±2.7 


9.13±0.15 


12.66±0.10 




LH-06 


33.7 ±3.7 


9.03 ±0.19 


12.46±0.15 




LH-03 


53.5±7.7 


8.66 ±0.20 


13.23±0.21 




L-1 


39.1 ±4.9 


8.79±0.20 


12.64 ±0.17 


2.16±0.18 


L-9 


41.3±6.4 


8.90±0.17 


12.82±0.24 


2.41 ±0.25 


L-11 


37.5 ±4.3 


8.85±0.18 


12.59 ±0.15 


1.83±0.16 


L-14 


40.4 ±7.3 


8.63 ±0.28 


12.61 ±0.22 


1.94 ±0.22 


L-15 


33.9±4.3 


8.81 ±0.22 


12.34 ±0.15 


1.99±0.16 


L-17 


43.3 ±8.9 


8.61 ±0.26 


12.75 ±0.29 


1.76 ±0.29 


L-20 


46.9 ±8.8 


8.63 ±0.26 


12.93 ±0.25 


2.48 ±0.26 


L-21 


36.0 ±5.1 


8.89 ±0.22 


12.52 ±0.19 


2. 17 ±0.20 


L-22 


33.4±3.5 


8.97±0.18 


12.41±0.15 


1.56±0.16 


L-23 


32.8 ±5.7 


8.95 ±0.24 


12.36±0.26 


1.80 ±0.27 


L-25 


36.3±5.3 


8.79 ±0.22 


12.48±0.18 


2.07 ±0.19 


L-27 


36.3 ±5.6 


8.77 ±0.25 


12.47±0.20 


2.30±0.21 


LH-01 


27.1 ±2.0 


9.35±0.13 


12.16±0.11 




ENl-01 


45.6 ±7.7 


8.66 ±0.24 


12.89±0.23 




ENl-02 


43.2 ±9.4 


8.65 ±0.29 


12.77 ±0.30 




ENl-04 


33.6 ±6.2 


8.89±0.28 


12.38±0.26 




EN2-01 


34.5 ±4.9 


8.82 ±0.23 


12.39±0.18 





Note. — All quantities were fitted using a single-temperature model with 
13 = 1.5 and assuming a 3kpc emission diameter. The dust masses assume 
k(850 /im)=0. 15 m^ kg"' . Uncertainties are la total errors of the fits to data, 
which do not include the uncertainties in the redshift values. 



Table 7 

Realistic Characterizations 





Tc 


logMrf 


logL 




ID 


(K) 


(Mq) 


(Lq) 




LH-02 


34.6 ± 1.9 


n A"7 1 A 1 T 

y.u/±o.i3 


1 T A~7 1 A AO 

lz.9/±U.U8 




LH-06 


33.5 ±2.9 


s.86±0.18 


12.93±0.12 




LH-03 


45.0 ±3.4 


5.71 ±U.12 


13.44 ±0.13 




L-1 


35.6±2.8 


oil 1 A 1 


13.04 ±0.12 


T 1 A 1 C 

1.15 ±U.15 


L-9 


34.8 ± j.O 


O Af 1 All 

o.yi ±u.ii 


12.9y±0.15 


O 1 A 1 


T 11 


33.z±z.4 


o.o4±U.l J 


1 o m _i_ n 1 1 

iz.y 1 it u. 1 i 


1 QZ _L All 


L-14 


33.5 ± 2.7 


o _i_ A 1 ^ 

o.oy ±U.lo 


1 T AT 1 n 1 o 

12.93 ±0.12 


1 A "7 1 A 1 C 

1.97 ±0.15 


T 1 C 

L-15 


29.8 ±2.3 


8. 78 ±0.16 


1 O T T 1 nil 

12. 73 ±0.1 1 


o no 1 n 1 y1 

2.08 ±0.14 


L-17 


34.1 ±2.9 


8.71 ±0.14 


12.96 ±0.14 


1.67±0.16 


L-20 


36.7 ±2.7 


8.77±0.13 


13.09 ±0.11 


2.34 ±0.14 


L-21 


31.1 ±2.4 


8.91±0.15 


12.80±0.12 


2.15±0.15 


L-22 


33.1 ±2.7 


8.79±0.21 


12.91±0.11 


1.76±0.15 


L-23 


33.0±2.7 


8.78±0.21 


12.90 ±0.14 


2.04 ±0.17 


L-25 


32.7 ±2.7 


8.74 ±0.17 


12.89±0.13 


2.17±0.16 


L-27 


31.9±2.7 


8.75 ±0.17 


12.84±0.13 


2.38±0.16 


LH-01 


27.7 ±1.7 


9.21 ±0.15 


12.58 ±0.09 




ENl-01 


38.5±3.1 


8.70 ±0.15 


13.17±0.13 




ENl-02 


34.7 ±3.0 


8.75±0.15 


12.99 ±0.14 




ENl-04 


32.0±3.0 


8.77 ±0.22 


12.85±0.15 




EN2-01 


32.5 ±2.9 


8.70±0.18 


12.88±0.13 





Note. — Same as Table|6] except that properties are derived for a power- 
law temperature distribution with 7 = 6.7, and a most-likely value of 2kpc 
is assumed for the emission diameter. We used K(850/jm)=0. 15 m- kg"' and 
(5=1.5 in deriving the quantities. 



SWIRE/MAMBO sources, and the classical SMGs. The re- 
sults from the single-T fits are statistically indistingu ishable 
from the values found for classical SMGs by Kovac Tet alj 
(|2006). As such, we tentatively conclude that their sam- 
ples consist predominantly of classical SMGs, although with 
a more restricted range of redshifts (z~2.05±0.30) than the 
purely submillimeter-selected samples. The narrow redshift 
distribution is a direct product of our selection criteria, which 
favor objects close to zss2. Since the radio co-selected SMG 
population peaks around z~2.3 (Chapman et al. 2003, 2005?), 
we conclude that two-pronged SWIRE selection (seeking the 
24 //m bright population of 5.8 /im peakers) is effective in 
identifying a representative subset of SMGs at mid-IR wave- 
lengths. 

5.1. Radio-FIR Correlation 

Our dust fits point to a tight correlation between IR and 
radio luminosities. This is witnessed by the consistent q val- 
ues in Tables |6] and |7] which were calculated according to 
Equation ( fT2] ). Such correlations are we ll known for various 
local galaxies (e.g. [H elou et al. 1988; Condon & BrodericM 
119911; ICondonlll992i: .Yun et al...2001.) . and extend to vari- 
ous populations into the distant universe (see [Ga rrett 2002 
lAppleton et alll2004t iBeelen et al.'l2006t iKovacs et al.ii2006l) . 
A flux density-base d measure o f the correlation was origi- 
nally formulated bv Helou et al.l ([1985) using 1.4 GHz radio 
and IRAS 60 /im and 100 /zm flux densities. However, since 
we do not have the data in the rest-frame IRAS bands to use 
with the original definition, we can either guess these from our 
models, or redefine the correlation in terms of the quantities 
we can measure more directly. 

The existence of the correlation is not surprising as both the 
synchrotron emission and the IR luminosities are driven by 
the short-lived population of high-mass (>8M0) stars. The 
synchrotron emission arises when relativistic electrons, which 
are produced by supernovae, decelerate in the magnetic fields 
of galaxies. At the same time, the IR heating is dominated by 



the luminosities of massive OB stars. With short-lived mas- 
sive stars responsible for both, measures based on luminosi- 
ties, not flux densities, are expected to capture the essence of 
the correlation more accurately. Unlike many prior studies 
of the correlation, we have the advantage of well-constrained, 
precise luminosity estimates. Accordingly, we can express the 
logarithmic correlation constant as 



? = log 



,4 GHz 



(12) 



with an appropriately chosen normalization frequency v„ 
for the IR luminosity. The rest-frame radio luminosity density 
^x.AGUz is calculated from the observed flux density S'i^ghz 
using: 



'Ci.4GHz = 47rD£ (1 +z)^ " 5i.4GHz 



(13) 



in terms of the luminosity distance Dl and radio spectral 
index a. The analysis of the SWIRE/MAMBO sample uses a 
values quoted in F09, while for the local starbursts we fitted 
a along with the other parameters. 

It is practical to choose the normalization frequency f„ 



2.6 
2.4 
2.2 
2 
1.8 
1.6 




Figure 4. Radio to (far-)infrared correlati on. We plot the q values for our 
sample, computed according to Equation 1121 using the luminosities from 
Table|7]and i/,i=9THz. Also shown ai'e the best-fit ensemble average value 
of 2.08±0.04 (solid line with shaded band), and the estimated 2a underlying 
scatter of data around it (dotted lines). 
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Table 8 

Comparison to Classical SMGs 









log Ml 


logMrf 


logLi 


logL 




Sample 


(K) 


(K) 




(Mq) 




a©) 


Tpk 7 


local stai'bursts" 




34.7±3.6 




7.33 ±0.62 




11.82±0.51 


2.42 ±0.35 0.3 ±0.2 7.22 ±0.09 



L09 36.1±5.8 34.4±4.2 8.96±0.18 8.85±0.12 12.59±0.28 12.97±0.22 ■■■ 1.4±0.5 

"¥& 36.5 33.1 8.83 8.80 12.54 12.91 2.08±0.16 1.1 ±0.2 

L09 + F09 35.5±2.2 33.5±2.1 8.90±0.07 8.82 12.56±0.05 12.93±0.12 ■■■ 1.2±0.4 6.71±0.11 

SMGs 34.1 ±9.5 (30.4±8.5)* 8.96±0.32 (8.85±0.18)* 12.60±0.48 (12.76±0.61)* 2.12±0.12 (1.2±0.5)* 

SMGs(1.5<z<2.5) 32.1 ±5.6 (30.8±5.5)'' 9.07±0.26 (8.87±0.11)* 12.57±0.30 (12.82 ±0.35)* 2.05±0.21 (1.2±0.4)'' 



Note. — A comparison of our 5piYz«r-selected samples to the local starbursts, and the classical SMGs of Kovacs et al. 1 2006). Indicated are the ensemble 
average values and 1 a intrinsic scatters for the characteristic dust temperature (i.e., single-T temperature) Ti, cold-component temperature Tc, single-T dust mass 
Ml , dust mass M^, single-r luminosity L\ , total (IR) luminosity L, the (F)IR-radio correlation measure q, and the characteristic opacity Tp^ near the peak of the 
emission. Only mean values are shown, when the intrinsic scatters are not detected. The last column shows the average mass-temperature distribution index 7 
and its rms uncertainty. The dust masses assume K(850/tm)=0. 15 m^ kg"' . 

" These are the local starbursts of Table[5] which provide the SWIRE starburst templates IBertall200^ : IPolletta et alJ2007l) . 
* Based on our own multi-r fits to the lKovacs et alj )200(i) data and assuming 7=6.7. 



such that the resulting q values can b e compared directly 
to the widely used original measure of lHelouet all CHI, 
which assumed typical temperatures of 40 K. We pick i/„ such 
that the definitions match exactly at that temperature, giv- 
ing i'„^>4.52THz used with Li (see iKovacs et al■ll2006^ . or 
I'm— >-9.0THz used with L assuming 7=7-2. Co nsequently, our 
q values closely mimic the lHelou et alJ (Il985h definition. The 
correlation measures q for our single- and multi-temperature 
fits are listed in Tables |6] and |7] respectively. 

A differ ent definition was suggested recently by 

llvison et"an (1201 Oh . who normalize IR luminosities, mea- 
sured in the 8-1000 /im band, with zy„=3.75 THz intended for 
the 40-120 ^m band by iHelou et"aLl (11985b . Their measure 
no longer offers a direct comparison to many earlier studies, 
which is why w e prefer our s . Ho wever, our q values can 
be compared to llvison et"an (120101) after adding 0.38 (i.e., 
log[9.0/3.75]). 

We find best-fit average q values of 2. 11 ±0.08 with the 
single-temperature dust models and 2.08±0.04 for the power- 
law temperature distributions. The implied intrinsic scatters, 
over the measurement uncertainties, are estimated at 0. 18 and 
0.16, respectively. Both our mean values and intrinsic scatters 
(Figure |4]i are significantly below those in local stud i es, an d 
comparable to the SMGs a nalyzed by i Kovacs et al.l (120061) . 
and to the ^ir~2.40 found bv llvison et a l. (2010) for their cor- 
relation measure. The result further supports the conclusion 
that our sample primarily consists of classical SMGs, without 
significant excess heating either by an AGN, or by low-mass 
stellar populations (which are likely partially responsible for 
the higher average q values observed in the local universe). 



6. CONCLUSIONS 

Our main conclusion is that the mid-IR selection criteria of 
the L09 and £09, picking the bright 24 /im sources that are 
also 5.8 /im peakers, is effective in identifying a significant 
fraction of the SMGs around z—2, thus providing a means of 
enhancing our understanding of the elusive SMCj population 
through studies at mid-lR and shorter wavelengths. The FIR 
characteristics (dust masses, temperature, and IR luminosi- 
ties) of the SWlRE-selected sources are essentially identical 
to those of the classical SMGs at the same redshift. Addi- 
tionally, the SWlRE-selected sample exhibits a correlation be- 
tween radio and IR luminosities that is indistinguishable from 
that seen for SMGs. Our other conclusions are as follows. 



1 . To provide realistic FIR characterizations for the sam- 
ple, we developed new models for the SEDs of galax- 
ies with power-law mass distributions of temperature 
components. We demonstrated that such models de- 
scribe local starbursts extremely well, under a tightly 
constrained dust emissivity index /3~1.5 and mass- 
temperature index 7~7.2. The 7 values are consistent 
with what we expect for heating dominated by star for- 
mation. 

2. We find that the most likely IR spectral continuum 
shapes of the distant SWIRE-selected starbursts dif- 
fer only slightly from the local examples, with best- 
fit 72±6.7±0.1. The lower 7 values are consistent with 
star formation under the high optical depths (rpk— 1.2) 
in our sample. 

3. Relying on the stacked IRS spectra of IL09I for estimat- 
ing the 24 /im continuum contribution, we place the 
typical diameter of starburst activity in these SMGs 
above 1.2kpc (with 95% confidence), and find a most 
likely value of 2 kpc. 

4. The emission scales, the typical 7, and the consistently 
low q values all confirm that star formation is the main 
power source in the galaxies of the sample. 

5. One of our targets shows a close clustering of three 
or four 350 /im components, each of which con- 
tributes comparably to the observable 350 /im flux 
density under the 15"-30" resolution of most current 
(sub)millimeter surveys. A commonality of unresolved 
SMG multiplets could bias our understanding of this 
population. We need further studies to explore the im- 
portance of close multiplets in current and future sub- 
millimeter surveys. 
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Frequency (THz) Frequency (THz) Frequency (THz) 



Figure 2b. Individual power-law temperature distribution SED fits to the local starbursts. The data used in our analysis are shown with hollow red squares and 
error bars. The sources with insufficient radio data assume a=0.75 for the radio spectral index. 
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Figure 2c. Individual power-law temperature distribution SED fits to our sources. The data are shown with hollow red squares and error bars. The 24 fim points 
reflect the expected continuum contribution to the in-band flux densities, according to Figure[3] with an uncertainty of 0.25 dex. The sources without radio data 
assume a=0.15 for the radio spectral index. 



